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Abstract-The grey Verhulst model has been widely applied in oil 
production forecasting. However, there are many defects of 
parameters estimation in traditional grey Verhulst model. To 
overcome the defects, an optimized grey Verhulst model based on 
BP neural network is proposed. Firstly, the Verhulst model is 
mapped to a BP neural network, the corresponding relationships 
between grey Verhulst model parameters and BP network 
weights is established. Then, the BP neural network is trained by 
means of BP algorithm, when the BP network convergences, the 
optimized weights can be extracted, and the optimized grey 
Verhulst model based on BP algorithm can be realized. By 
applying this optimized grey Verhulst model to the national oil 
production forecasting, the experiment results show that the 
method is feasible and effective, especially, when the oil 
production increases according to the curve with S-type, not only 
higher forecasting accuracy can be obtained, but also the 
superiority and the features of grey system model can be reserved. 
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I. INTRODUCTION 

With the fast economic development, the energy demand is 
increasing year after year, especially the oil demand. The 
supply of oil is a hot problem. The condition of oil reserves in 
china is not good. Most oil field has been at the stage of 
mining metaphase or mining anaphase. The increase of oil 
production is very difficult. Facing the rising oil demand 
increasingly, the forecast of oil production scientifically and 
exactly in the future has realism significance to the plan of oil 
demand and oil stratagem repertory [1, 2]. There are many 
methods for oil production forecasting, and the grey system 
model is one of them The GM (1, 1) model is widely 
discussed and studied in the grey system. The purpose of GM 
(1, 1) model is to work on system forecasting with poor, 
incomplete or uncertain messages. In the prediction of oil 
production, we usually use the grey GM (1, 1) model. 
However, it is imperfect when the oil production increases in 
the curve with S-type, or the increment of oil production is in 
the saturation stage. In this case, the oil production forecasting 
error of grey system model will become larger and the result is 
unaccepted in the real world. The grey Verhulst model is a 
special kind of model within the grey system. The grey 
Verhulst model [2-9] is used to forecast the raw sequence 
growth in S-type or growth being saturated. Therefore, in 
recent years this model has been widely applied in oil 
production forecasting, especially, when the oil production 
increases according to the curve with S-type [1, 2, 3]. 
However, there are many defects of parameters estimation in 
traditional grey Verhulst model [2-8].First, the traditional 
model need to face the reasonable selection of background 
value. On the other hand, in order to obtain forecasting model, 



in these improved models, parameters which obtained from 
the grey difference equation is substituted into the grey 
differential equation, thus we also can not avoid the jumping 
errors from the difference equation to differential equation in 
traditional grey modeling [2, 10-12]. This parameter 
estimation method which obtains through the indirect method 
may not guarantee that the value and the actual value of the 
gray prediction error are minimum. Generally, from the 
evaluation criteria of the model fitting, when solving grey 
Verhulst model model, the model should directly take 
minimizing the error of the predicted value and the actual 
value as the criterion. 

To consider the forecasting of the raw sequence growth in 
S-type or growth being saturated, from the data fitting's 
viewpoint, this study uses BP neural network model to 
overcome the defects of grey Verhulst model. The grey 
system's primitive discrete data is fitted by a continuous model 
[13] which has the same form with grey Verhulst model's time 
response type, the model parameters are trained and optimised 
by means of BP neural network, and then the optimization grey 
Verhulst model is realized. 

The remainder of this paper is organized as follows. In 
Section 2, the grey Verhulst model is introduced. In Section 3, 
the Verhulst model is mapped to a BP neural network, the 
relationships between rey Verhulst model parameters and BP 
network weights is established, where BP network algorithm 
is used to optimize grey Verhulst model parameters. The 
experimental results and discussions are presented in Section 4. 
Finally, conclusion is drawn in Section 5. 

II. GREY VERHULST MODEL 

Grey Verhulst model constructing process is described 
below: 

Denote the nonnegative original sequence by 

X (0) =(x (0) (l),x (0) (2),...,x (0) («)). (1) 

The accumulated generating operation (AGO) [1] 
formation of X <0) is defined as : 

X m =(x m (l),x m (2),--;x m (n)), (2) 

k 

where x m (l) = x (0) (l),x m (k) = X* (0) (;), it = 2,3,-, n. 

7=1 

The grey Verhulst model can be established by 
constructing a first order differential equation for x a > 
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dx (1) 
dt 



+ ax m =b(x m ) 2 



(3) 



The solution of Eq.(3)can be obtained by using the least 
square method, that is: 



x m {k + l)- 



ax (0) (l) 



bx m (1) + {a - bx m (l))e" 



(4) 



where 













[a,b] = (B T BY 1 B T Y and B = 




(^) 2 


Y = 






-z m 






x (0) 

"Si 



z t m = Q5{x w (k) + x m (k~\))- 

Applying the inverse accumulated generating operation 
(IAGO) [1], the grey Verhulst model fitted and predicted 
values are obtained by 



x m (k) = x m (k+Y)- x m (k) (k = 2,3,-- •) 
x (0 \l) = x (0 \l). 



(5) 



Assume that the original sequence X i0) itself increases in 
the curve with S-type or the increment of the original sequence 
is in the saturation stage, the original sequence can be taken 

as » the IAGO of the original sequence can be taken 
asZ (0) . 

HI. GREY VERHULST MODEL BASED ON BP NEURAL NETWORK 

Assume that X (<)) increases in the curve with S-type or 
the increment of the original sequence is in the saturation stage, 
to consider the following non-linear function, 



where a,b,x (0) (Y)^ undetermined parameters. 



bx m (l) + (a-bx i0 \iy)e" ( - 



(6) 



Compared with grey Verhulst model, the function value of 
formula (6)at t ( t = 2,3,---,n ) can be taken as the 

predicted value in the corresponding time t of the grey 
Verhulst model. In order to obtain a better solution, this initial 
value does not take the first data of the original sequence data, 
but as a parameter to be determined by data fitting. Hence, the 

original discrete data sequence X ( °'can be fitted by formula 
(6). The key question is how to estimate the parameters by 

original discrete data X (0) . 

Because the error function is nonlinear, we now resolve the 
above problem by means of BP network which can 
approximate nonlinear function with any precision [14]. 

First, the formula (6) is transformed as follows: 

1 



y{t) = 



l 



1 ^v a( ' _1) 

a V>(1) a 



(V) 



(8) 



y(t) a x (u, (l) a 



Since e " (t l) > 0, therefore, 



1 b e MH) b 1 
- + — 



y(t) \ x (1) a (l + e ) a \ + e J 



1 



V x (1) a 
1 b 



1 



+ 

a\ + e 



1 126 1 

+ 



\ x (l) a x (1) (l + e ) a i + e J 



+ e ) (9) 



The above equation is mapped to the BP neural network, 
the network structure is described as Fig. 1: 




LA LB LC LD 



Figurel Grey BP neural network structure 

The corresponding relationships between the network 
weight and the grey Verhulst model parameters are established 
as follows: 

1 



<^n - a > a> 7 , = — - 

x (0) (l)' 



2b 



6) 3l =a> 32 =l + e 



-0(1-/) 



(10) 



The threshold value of LD level is 

And the LB level neuron's transfer function is taken as 

1 

-, there is a high-gain 



Sigmoid type functions f(x) = 



l + e~ 



area which ensure that the network eventually reach a steady 
[13] .Other neuron's transfer function is taken as an e linear 
function f(x) = X ■ 

Through the above method, the grey Verhulst model is 
mapped to a BP neural network; figure 1 shows this neural 
network structure for grey Verhulst model is simple. The 
corresponding relation between the grey Verhulst model 
parameters and BP network weights is established, then the 
grey Verhulst model parameters estimation problem is 
transformed into the weights of neural network optimization 
problem. The BP neural network is trained by means of data 

1 , when the BP network 



sets 



(l-t, 



* <0, (0 



),t=l, 2,---,n 



convergence, the optimized model parameters can be obtained 
[13], and the optimization modeling for the grey Verhulst 
model based on BP algorithm can be realized. The model can 
avoid the jumping errors from the difference equation to 
differential equation, and overcome the defects of the grey 
Verhulst model. 
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IV. EXPERIMENT RESULTS AND DISCUSSIONS 

A. Experiment Design 

In this section, to validate the performance of the grey 
Verhulst model based on BP neural network. The optimized 
grey Verhulst is applied for national oil production forecasting, 
the primary data is shown in Table I [1]. 

TABLE L National oil production from 2002 to 2009 [Unit: 
million tons] 



Year 


200 
2 


200 

3 


200 
4 


200 
5 


200 
6 


200 
7 


200 
8 


200 
9 


Serial 
number 





1 


2 


3 


4 


5 


6 


7 


Product 
ion 


16.7 



16.9 
5 


17.5 

8 


18.1 

3 


18.4 
7 


18.6 

3 


19.0 



18.9 
4 



The actual value data presents the saturation characteristic; 
therefore, it is suitable for grey Verhulst modeling. According 
to the proposed method in this paper, the BP neural network 
with three hidden layer is established. The network structure is 
described as Fig.l. The first level neuron's transfer function is 
taken Log-Sigmoid type function; other neuron's transfer 
function is taken as a linear function. The maximum training 
epoch is M=300 , the allowance permissible error is 
E = 0.00000005 ■ The LM (Levenberg -Marquardt) is taken as 
the training algorithm; the learning rate is dynamically 
determined by LM algorithm. To use the Matlab6.5 
programming, when the network achieves the accuracy 
requirement, we can obtain the optimized grey Verhulst BP 
neural network model. 



Performance is 4.56462e-008, Goal is 5e-008 




2 3 4 
6 Epochs 

Figure2 Training errors 

Fig. 2 shows the training times of this method only use 6 
steps to reach the minimum error. However, the traditional BP 
network not only need a large amount of data to train, but the 
theoretical guidance of network settings is also lack moreover 
the training time is much greater than this article method [14], 
Many experiments show that the network convergence epochs 
of the method is generally about 10. Experiment shows the 
method combines the characteristics of the small sample and 
poor information of the Verhulst model, and displays the 
characteristics of strong nonlinear approximation and fault- 
tolerant capability of the neural network. 

B. Empirical Results and Error Analysis 



Compared with the traditional grey Verhulst model and other 
method in the paper [1], the fitted value of different models and 
the actual value are shown in Fig. 3 for comparison. 




Serial number 

Figure3 Production of crude oil from 2002 to 2009 fitted by different models 

Mean absolute percentage error (MAPE) approach [15] has 
been recommended to validation. Mean absolute percentage 
error is defined as 



MAPE 



i i 



y, - y, 



(ii) 



where y. is the actual value, and y\ is the predicted value, 



y, - y t 



y, 



is absolute percentage error of y ; . The Experiment 



results are shown in Tab. II . 
TABLE E. 



The forecasting relative errors based on 
different models 



Model 


Grey Verhulst 
model 


Paper[l] 


Ihis paper 


MAPE % 


0.76 


0.75 


0.32 



According to results listed in Tab. II , compared with the 
traditional model, the precision in this article enhances57.89%. 
Compared with other method in the paper [1], the precision in 
this article enhances 57.33%. Analyzing the errors of 
competitive models by means of comparing the MAPE, we 
conclude that the method in this paper is a better fitted 
prediction model. Furthermore, according the criteria purposed 
by Lewis [15]. The model is better than the other forecasting 
models because of lower MAPE. 

V. CONCLUSION 

1. From the data fitting's viewpoint, the grey Verhulst 
model can be optimized by means of BP algorithm. 

2. The optimized grey Verhulst model in this paper can 
overcome the defects of the grey Verhulst. 

3. The experiment results show that the optimized grey 
Verhulst for oil production forecasting is applicable, especially, 
when the oil production increases according to the curve with 
S-type, not only higher forecasting accuracy can be obtained, 
but also the superiority and the features of grey system model 
can be reserved. 
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